hcp-ph/980529l| SINP/TNP/98/11 

Gravitational coupling of neutrinos in a medium 



Jose F. Nieves 

Department of Physics, P.O. Box 23343 
University of Puerto Rico, Rio Piedras 
Puerto Rico 00931-3343 



p. 

> 

X 



Abstract 



Palash B. Pal 

^ ! Saha Institute of Nuclear Physics, 1 / AF Bidhan-Nagar 

; Calcutta 700064, India 

May 1998 

< 

In a medium that contains electrons but not the other charged leptons, such as normal matter, the 
■ gravitational interactions of neutrinos are not the same for all the neutrino flavors. We calculate the 

^ ' leading order matter-induced corrections to the neutrino gravitational interactions in such a medium and 

consider some of their physical implications. 

a^ 
<^ 

^ ! 1 Introduction 
O 
oo 



When neutrinos propagate through a medium, the effects of the background particles can influence their 
properties in important ways. The possible consequences of these effects have been the subject of continuous 
research activity in recent years, largely motivated by their implications in a variety of astrophysical and 
cosmological contexts. Some examples of the latter include the original MSW mechanism ^ and its variations 
p> il in the context of the Solar Neutrino Problem |Q , and the explanation of the large birth velocities of pulsars 

^ ■ in terms of the asymmetric emission of neutrinos from the cooling protostar, which is produced by the 

matter-enhanced neutrino oscillations biased by the supernova's magnetic field |^ . 

Among the various approaches that exist to the study of the matter effects on the propagation of neutrinos, 
the methods of Quantum Statistical Field Theory (QSFT) have proven to be very useful ones. These methods 
5^ ' have been applied to reproduce the Wolfenstein formula for the index of refraction of neutrinos in matter 

I and to obtain corrections to it 0, ||, |6[ 0] . In addition, they have been used to determine the electromagnetic 

properties of a neutrino that propagates in matter by means of the one-loop calculation of the electromagnetic 
vertex function induced by the neutrino interactions with the background particles Furthermore, it was 
observed in Ref. [|j that, in the presence of a static magnetic field, the effective electromagnetic interactions 
of the neutrinos produce an additional contribution to the neutrino index of refraction which modifies the 
condition for resonant oscillations in matter. This effect is in fact the origin of the possible explanation of 
the large birth velocities of pulsars mentioned above. 

In all these situations, a common theme has been the observation that the universality of the neutrino 
interactions is broken due to the fact that normal matter contains electrons but not the other charged leptons. 
Therefore, in a medium such as the Sun or a supernova, the electron neutrinos on one hand, and the tau 
and muon neutrinos on the other, are affected in different ways. This fact implies that in a such medium 
the GIM mechanism is not operative, and it explains why the radiative decay of a neutrino in a medium is 
greatly enhanced compared to the corresponding rate in the vacuum pO{ . 

The observation that the universality of the neutrino interactions is broken by the background medium is 
a special case of a more general concept. It is the notion that while the fundamental Lagrangian of the theory 
is invariant under certain symmetry operations, a given background medium may not be invariant under 
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the same operations. Thus, for example, normal matter is neither CP nor CPT asymmetric, a fact that 
has interesting consequences for the propagation of photons in a (chiral) medium such as a gas of neutrinos 
pl| . Taking this notion to the extreme, we can even think of the medium as breaking Lorentz invariance 
since it specifies a preferred frame of reference; i.e., that in which the medium is at rest. In the context of 
QSFT, this breakdown of Lorentz invariance is traded by an explicit dependence of the calculated physical 
properties of the particles on the velocity four-vector of the background medium. In this way, for example, 
the photon acquires a gauge invariant mass, a left-handed neutrino acquires a chiral-invariant mass and a 
Majorana neutrino acquires electromagnetic dipole moments. 

The main observation of the present paper is that, in the presence of medium, the breakdown of the 
universality of the neutrino interactions includes their gravitational interactions as well. In analogy with the 
fact that, in a medium that contains electrons but no muons or taus, the electron neutrinos have different 
electromagnetic interactions than the muon or tau neutrinos, their gravitational interactions also differ. 
In this work we determine the effective gravitational interactions of neutrinos in a matter background, by 
calculating the one-loop contribution to the neutrino stress-energy tensor, which is the gravitational analog 
of the electromagnetic current. 

We wish to clarify at the outset the following point. The possibility that the observed deficit of solar 
neutrinos can be ascribed to neutrino oscillations that are driven by flavor changing interactions of grav- 
itational origin has been proposed previously In these works, it is postulated that the gravitational 
interaction of neutrinos has the same form as the standard one, but with different coupling strengths for 
each neutrino flavor. This assumption violates the equivalence principle in a fundamental way, and in fact 
destroys the symmetry that makes it consistent to assume the existence of a massless graviton in the vacuum. 
On the other hand, no such fundamental breaking of the equivalence principle is assumed in the effect that 
we describe in the present work. Our calculations are based on the standard model of particle interactions, 
together with the commonly accepted linearized gravitational coupling of fermions via the stress-energy ten- 
sor, with a universal coupling strength. The non-universal character of the effective gravitational coupling 
of the neutrinos, whose calculation is the aim of this work, is a consequence of the flavor asymmetry of the 
medium and not of the fundamental Lagrangian, as we have emphasized above. 

We have found convenient to organize the presentation in the following manner. In Section 2, we outline 
the linearized theory of the gravitational coupling of fermions, including the interaction terms with the W 
and Z gauge bosons. This has been necessary because the formulas that are commonly used and/or quoted in 
the literature are valid when the fermions that couple to the graviton are on their mass-shell. Therefore, it is 
not appropriate to use them for the internal fermion lines in the one-loop diagrams that we need to consider. 
Armed with the preliminaries of Section |^, the one-loop diagrams for the induced gravitational interactions 
of the neutrinos in a background composed of electrons and nucleons are calculated in Section ||. There we 
obtain the effective gravitational vertex function as an integral over the momentum distribution functions 
of the background particles. We conclude that section with some general remarks about the consistency 
of the one- loop formulas, and in particular we show explicitly that our result for the effective gravitational 
vertex of the neutrinos, including the background-induced terms, is transverse as it should be. Using the 
one-loop formulas for the vertex function as a starting point, we determine in Section |^ the modification to 
the neutrino index of refraction in a medium in the presence of a static gravitational field. There we consider 
also some possible applications of these results, giving some of the details in the Appendix, and Section ^ 
contains our conclusions. 



2 Tree-level gravitational couplings 

It is well known that the linear, or weak field, approximation to the metric theory of gravity is equivalent 
to a quantum particle description of gravitational interactions in which the graviton emerges as a spin-2 
quantum field coupled to the stress-energy tensor. The formulas for the gravitational vertices is given and 
discussed in many textbooks, at least for the common cases of spin 0,1/2, and spin-1 (i.e., photons) particles. 
However, the formulas that are customarily quoted are given for the case in which the particles are on their 
mass-shell, which is sufficient for the applications that have been considered in the literature, but not for 
the present one. This is particularly true for the fermions. 

In the one-loop calculation that we are considering, there are diagrams in which the graviton couples 
to internal electron lines in a loop. For them, the on-shell form for the gravitational vertex is not valid. 
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Furthermore, there are other interaction vertices that are unique to our case in hand. In this section we 
consider all the couplings that relevant to our calculation in detail. We define the action A in presence of 
gravitation by 



J d^xC, 

and then look at the different terms in C which will be relevant for us. 



(2.1) 



2.1 Fermion couplings 

The Lagrangian for a free Dirac fermion of mass m f can be written in the explicitly Hermitian form 



(/) 



1 

--07'' 9^7/1 + h.c. 



(2.2) 



The Lagrangian in the presence of gravitational interactions is obtained from this by making the replacement 
l^d^ ^"''Da, where the 7" denote the ordinary gamma matrices and 



(2.3) 



is the gravity-covariant derivative, with a 
defined in a way that 



ab i r_,a 



[7°i7 ]i and Va'^ being the vierbein vector fields. These are 



gti.iyVa^Vb" = r^ab , 



(2.4) 



where g^i, is the space-time metric and r]ab is the flat space metric, which we take to be diag(l, —1, —1, — 1). 
The spin connection coefficients are given by 



(2.5) 



where ^x^i^ are the Christoffel symbols. This formula can be expressed explicitly in terms of the vierbein 
fields in the form 

Wah7 = ]^{va^{^^Vb^,- d^Vb-y) + Vb'^id^.Va^ ~ d^Va^,) Va''vb''vc^{d^v'' ^- d^v"^)} . (2.6) 

In addition to these changes, we have to include the determinant of the matrix Uq^, which we denote by 
det(w), as an overall factor. 

Thus, in the absence of all interactions except the gravitational ones, the Dirac Lagrangian is given by 



= det(w) 



i^l'^Va'' - -UJbcy^ + h.c 



(2.7) 



With this construction, the term in the action corresponding to Cg is invariant under general coordinate 
transformations, labeled by the greek indices, and also under local Lorentz transformations labeled by the 
latin indices. 

The linear theory of the gravitational couplings is obtained by assuming that, for weak gravitational 
fields, we can write 



(2.8) 



and treat the second term on the right side as a perturbation. The quantity k is related to the Newton's 
constant G through the equation 



(2.9) 



in order that identified with the graviton quantum field, has the correctly normalized kinetic energy term 
in the Lagrangian. The vierbeins cannot be determined uniquely from Eq. (2.4). They can be determined 
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only up to a local Lorentz transformation which, however, would leave the action invariant. Therefore, we 
can choose in particular 



which in turn gives 

Note that this also implies that 



det(w) = 1 + Kr]^^h'''' 



since the matrix is the inverse of the metric g^j^ and from Eq. (2.4), 



(2.10) 
(2.11) 
(2.12) 
(2.13) 



We now substitute these relations into Eq. (2/7), and keep only the terms that are at most linear in k. Once 
this approximation is made, it is no longer necessary to distinguish between greek and latin indices (since 
any difference between them would contribute only to hig her order in k) and therefore we write everything 
in terms of the greek indices from now on. From Eq. (2_^) we obtain 



UJn 



2 {di/hf^ij 



(2.14) 



and it then follows that the term involving the cr-matrices in Eq. (2.7) drops out because it is proportional 
to the quantity 

{J^,ax,} + {^I^ly) ^0. (2.15) 
From the other terms we then obtain the gravitational interaction Lagrangian of the fermion in the form 



C 



Uf) 



where the stress-energy tensor operator T^{.^ for the fermion field is given by 



-tlj{x) Yf^dy + "fyd^] il){x) + h.c. 



Vfii'^O ^ (^) ■ 



(2.16) 



(2.17) 



From Eqs. (2.16) and (b.l7) it follows that the term correspondmg to the gravitational fermion vertex m a 

If) 

Feynman diagram is —inVi^i , where 

1 



(2.18) 



In the above considerations we have assumed that the fermion / is a Dirac particle with a given mass. 
On the other hand, it is easy to infer by inspection that similar arguments yield the formula 



V^-^Hk, k') = \ [7p(fc + k').+ 7.(fc + k')^\ L - ir,^. + ^']L 



(2.19) 



for the case of a chiral, left-handed neutrino. 



Notice that if we take the matrix element of the operator tI£\o) between on-shell fermion states, with 
incoming and outgoing momenta p and p' respectively, the term proportional to 77^^ in Eq. (2.17) gives no 
contribution and we obtain 

fip') TyXo) f{p))^]uip')H.{p + P%+lAP + p')^}u{p), 



(2.20) 

which is the expression that is quoted in textbooks p4[ |. However, for the purposes of the one- loop calculation 
that we carry out in Section |^, we need to use the vertex for off-shell fermions given in Eq. ( 2.18| ). In 
particular, as we will show in Section ^ the one-loop calculation of the effective gravitational vertex of the 
neutrino gives a result that satisfies the transversality condition provided that the term proportional to Ty^jy 
in Eq. ( ^.1^ ) is included in the calculation of the loop diagrams. 
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2.2 Fermion and W boson couplings 

The interactions that drive the effective gravitational vertex of the neutrinos, are the standard weak inter- 
actions with the particles of the background. Let us consider the charged-current interactions first. In the 
presence of a gravitational field, the interaction term 6^7'^ W^j^l in the Lagrangian is modified according to 



= det(w) 



v2 



In the linear approximation given by Eqs. ( |2.1C| ) and ( 2.11 ), this becomes 



(cc) 



,(cc) 



where C^q'^'^ is the standard charged-current interaction Lagrangian 



V2 



(2.21) 



(2.22) 



(2.23) 



while 



C 



(cc) 



V2 



V^iuWaCLl'^UL + - {W^^eL^luVL + Wi^ClIp^Vl) + h.c. 



where 



-K-^ a^^Ap eL-l^UL WPh^"" + h.c. , 



(2.24) 



(2.25) 



2.3 W boson couplings 

The gravitational vertex for photons is well known. For the W boson, the only difference comes from the 
mass term in the kinetic energy part of the Lagrangian. Indeed, from the usual expression for the kinetic 
energy of the W vector boson, 



£w = det{v) j-ii^;,!^^'' - M^w;w^^^ , 

and then making the substitutions given in Eqs. (U) and ( pd| ), we arrive at 

£(W) ^ £(W) 

where £('^) is the standard form of the kinetic energy for the free W boson, while 



{WW) 



W*^pW^P - M^W*W^ rj^, + {W;j¥^ - M^W;W, + h.c.) 



(2.26) 



(2.27) 



(2.28) 



gives the gravitational vertex. 

Of particular importance for us are the terms proportional to the W mass. The reason is the following. 
In the analogous calculation of the charged-current contribution to the electromagnetic neutrino vertex in a 
medium, the dominant term is of order 1/M^ and it arises from the diagram that corresponds to diagram 
(A) in Fig. ^, in which the photon couples to the electron line in the loop. The diagram in which the photon 
couples to the internal W line, corresponding to diagram (B), is of order 1/M^, and therefore it is negligible. 
In the present case, because of the presence of the term that is proportional to the W mass in the coupling 
of the W to the graviton, diagram (B) gives a contribution also of order that must be taken into 

account. As far as these terms only are concerned, we can then replace Eq. (2.28) by 



(WW) 



(2.29) 



where 



(2.30) 



5 




2.4 Including the neutral-current couplings 

The charged current interaction discussed above are relevant only for the i^e's. The v^^s and the Vr^s will 
interact with the electrons only through the neutral current. Moreover, a normal background contains 
nucleons as well, which interact with all neutrinos via the neutral current. To take these interactions into 
account, we consider the neutral current couplings of the Z-boson. In absence of gravitation, these are 



4"^) = -gzZ^ 



—e.fijT f—e,p,n 



(2.31) 



where, in the standard model. 



gz 

-Xp = Xn 



g/{2 cos 9w) 

1 9 

2 sin 9w 



Xn 
Ye 



^9A. 



(2.32) 



with gA ~ 1.26 being the rcnormalization constant of the axial-vector current of the nucleon. It is also useful 
to remember that, in the standard model, 



2 2 

9z _ 9 _ jT^i^ 



(2.33) 
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Figure 2: Z-exchange diagrams for the one-loop contribution to the gravitational vertex of any neutrino flavor 
[i — e,iJ,,T) in a background of electrons and nucleons. 



Following the arguments that were used for the charged current, we find that in presence of gravitation, 
these interactions will be modified to 



/•(nc) _ /-(nc) ^(nc) 
'-g — '-0 '-h ' 

where £q"'^'' has been given above in Eq. ( ^.31 ), and 



C 



(nc) 



-ngz afj_^\p 



—e.fi,T /— e,p,n 



(2.34) 



(2.35) 



where a^j^^xp has been defined in Eq. (2.25). Also, there are terms which comes from the Lagrangian of the 
pure Z-boson which give interactions with the graviton and are proportional to Mf . These are 



(2.36) 



where a'^jjXp is the tensor defined in Eq. (2.30) 



3 Induced gravitational vertex of neutrinos 

We now consider the background contributions to the neutrino gravitational vertex. We denote the proper 
vertex function for off-shell neutrinos, including the background-induced contributions, by T^i^} {k,k'). It is 
is defined such that the matrix element of the total stress-energy tensor operator T^j^ (x) , between incoming 
and outgoing neutrino states with momenta k and k' respectively, is given by 

(z.(fc')|f^.(0)|i^(fc)) =IZi(fc')r|rJ(fc,A:')«L(fc). (3.1) 

It is useful to divide the relevant one-loop diagrams into the two sets given in Figs. ^ and according to 
whether they involve the W or the Z boson. In the diagrams, and in the formulas that follow, q stands for 
the outgoing graviton momentum 

q = k-k'. (3.2) 

The additional diagrams in which the graviton line comes out from one of the external neutrino legs are 
not shown, since they are 1-particle reducible and do not contribute to r|it''. The proper way to take them 
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into account in the calculation of the amplitude for any given process, is by choosing the external neutrino 
spinor WL(fc) to be the (properly normalized) solution of the effective Dirac equation for the propagating 
neutrino mode in the medium, instead of the spinor represent ing the free-particle solution of the equation in 
the vacuum. This will be discussed in more detail in Section |3.4 

As commented earlier, the VF-exchange diagrams of Fig. |l] contribute only to the gravitational coupling of 
i^e. In contrast, the diagrams shown in Fig. |^ contribute equally to the gravitational vertex of all the neutrino 
flavors Ve.ii.T, and therefore are not relevant in phenomena that involve transitions between the standard, 
weak SU(2)-doublet neutrinos. However, in processes in which the so-called sterile neutrinos participate, 
these flavor-diagonal contributions are important. Moreover, we consider them first since the results for the 
Ty-diagrams can be easily obtained from the corresponding ones for the .Z-diagrams by making some simple 
substitutions. 



3.1 The Z-mediated diagrams 
3.1.1 Diagram (A) 

This contribution, which will be marked by the superscript 2A, is given by 



Ml 



E 

f=e,p,n 



(27r) 



Tr 



7A {Xf + y/75) iS^J'> (p'Wy,^ (p, p')*4^' (P) , (3.3) 



or equivalently 



f=e,p,n 



Tr 



7a(A/ + 17.75) ^si^\p'MP{p,p') iSi^\p) 



In these expressions we have introduced the shorthand 

p' = p-q, 



(3.4) 



(3.5) 



which will be used throughout the calculations, and we have neglected the momentum dependence of the Z 
propagator in anticipation of the fact that we are interested only in the terms that are linear in the Fermi 



constant Gp- Further, Sp is the thermal fermion propagator which is given by 



zS'-J\p) 



p^ — mj: 



where 



77/ (p) 



0{p ■ v) 



2tt5(p^ - m))'qf{jp) 



-p-v) 



g0{p-v-fj.f) I ^-I3{p-v+fj,f) I 



(3.6) 



(3.7) 



with (3 being the inverse temperature, /ij the chemical potential and u'' the velocity four- vector of the 



(3.8) 



medium |15| . In the frame in which the medium is at rest, v'^ has components 

= (1,0) . 

( f) ( f) m K A 1 

When the expressions for Sy!'{p) and S)^'ip') given Eq. {pM are substituted in Eq. (|3.4|), three kinds o:' 
term are generated, according to whether they contain none, one or two factors of 7]f. We are interested 
in the terms with at least one factor of 7y/, since that is where the background dependence of the neutrino 
gravitational vertex will come from. However, the terms with two factors of 77^ are relevant only in the 
calculation of the absorptive part of the amplitude, which we are not considering here, and therefore the 
important terms for us are the ones having a single factor oi rjf. Denoting their contribution to the vertex 

function by T'^j^^^ we then obtain 



ni"^ - ^^'^ E /7S3Tr[7A(A,+l>75)A,.] 




(3.9) 
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where 



Using for V^l the expression given in Eq. ( 2.18D , we have 



(3.10) 



(3.11) 



where 



a-ii' = (|5' + m/)-{7^(p + p')i'+7>^(P + ]3')AJ(l^ + ™/)- 



(3.12) 



When Eq. ( 3.11 ) is substituted in Eq. (3^), the terms that involve the factor of 7/^,^ are reduced very simply, 
and the result can be expressed in the form 



(3.13) 



where a\^^ is given by the same expression given in Eq. ( |3.9| ) but with the replacement A^j^ a^^, and 



with 



bf = aV2Gf 



6{p'^ — rn?f)rif{p) p ■ v . 



In arriving at Eq. (3.13) we have used the formula 

d^p 



5{p - mj)r]f{p)p^ = hfv^ , 



(3.14) 



(3.15) 



(3.16) 



and the relation V^Gp = g^/M^- 

It is convenient to introduce the particle and antiparticle momentum distribution functions 

where the upper and the lower signs hold for the particle and the antiparticle respectively, and 



The corresponding total number densities are given by 



and in terms of them Eq. ( 3.15| ) can then be written in the form 

bf = V2GF{nf +nj) . 



(3.17) 



(3.18) 



(3.19) 



(3.20) 



The expression for AJ;^ can be simplified by taking the traces and making change of integration variable 
p — > p + q in the term containing the factor of ?7/(p')- this way we obtain 



A(2) 



9z f^,\ 
Ml 



(2^) 

'XfN^^J,{p,q)-YfN;^J,{p,q) XfNl'J,{p,-q)+YfN^^J,{p,-q) 



2p-q 



q^ + 2p- q 



(3.21) 
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where 

It is useful to observe that these quantities satisfy the relations 



KJx'i-P'-<l) = -K^^(P^l)' 



(3.22) 



(3.23) 



from which other similar relations can be obtained. Using them, after performing the p° integration in Eq. 
( 3.21 ) we obtain finally 



^2 -2pf -q 



+ {q^ -q) 



Yf{ff + fj) 



gf2 -2pf -q 



(q^-q) 



(3.24) 



Further reduction of Eq. (B.24) is not possible without making some assumption about the conditions of the 
electron gas, in order to be able to carry out the required integrations. Before turning to that, we consider 
the contributions from the other diagrams, that must be added to Eq. ( 3.13| ) to yield the full one- loop result 
for the induced vertex function. 



3.1.2 Diagram (B) 

Because this diagram contains two Z-boson propagators, we retain only the terms proportional to from 
the ZZh vertex given in Eq. ( 2.36| ), since the rest will yield results proportional to For the same 

reason, we also neglect the momentum dependence of the two Z propagators in the diagram. In this way we 
then obtain 



Ml 



'Y{Xf + Yn^)iS^P{p)\ (3.25) 



for the contribution to the vertex function. For the background-dependent part, this yields 



(3.26) 



/ 



which, by taking the trace and using Eqs. (|3ll ) and ( [sl^ ), can be reduced to 

p/(2B) _ -h^Z) I ^T^.P 

= -b^^'>[i]f,^]!! + jf,v^+j^Vf,]L. 

3.1.3 Diagrams (C) and (D) 



(3.27) 



For these diagrams we use the vertices given in Eqs. (2.31) and (2.35). The 0(1/Mf ) term is given by 



(27r) 



rTr 



r{Xf + Yfj5)^si^\p) 



(3.28) 



(2D) (20 

and, since a^^Xp is symmetric in the indices A and p, it follows that Vin, — Vin, . Following the same steps 
that led to Eq. ( |3.27 ) we find for the background dependent part 



(3.29) 
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3.2 VF-exchange diagrams 
3.2.1 Diagram (A) 

We now consider the diagrams shown in Fig. |l|. We begin with diagram A, which corresponds to the 



expression 



*5 \ / ■ N 



(2^ 



where we have neglected the momentum dependence of the W propagator. Using the identity 

(7^i)M(7AL) = -(7^L)Tr(M7Ai) , 
which is vahd for any 4x4 matrix M, Eq. ( 3.30| ) can be written as 



w 



(27r)4 



Tr 



'tS^;\p')V^.{p,pysP{ph^L 



(3.30) 



(3.31) 



(3.32) 



This expression coincides with what is obtained from Eq. (3.4) by discarding the nucleon terms and then 
making there the substitutions 



g| ^ 9^ Y Y 



(3.33) 



With this observation, the result for the background-induced part of t']}^'^ can be immediately deduced to 
be 



IJ.V } 



(3.34) 



where 



[fe - /-) 



2£;e(27r)3 

+ (/e + /-) 



-2pe- q 



+ iq^ -q) 



<Abe,g) 

g2 ~2pe- q 



(q^-q) 



(3.35) 



3.2.2 Other diagrams 

The steps to follow in the calculation of the other diagrams of Fig. |l| are very similar to those for the 
corresponding diagrams of Fig. ||. Omitting the details here, the background-dependent terms are given by 
the formulas 



V2 



(27r)3 



V ) i-i^^a^.o^p) I -^Sip^ - ml)7^,ip)rL{^ + m,)^^L , (3.36) 
w 



(27r)3 



with E'li^"* = V'^^J . Using Eq. ( ^.16| ), after some straightforward algebra these reduce to 



/(C) 



p/(lB) 



p/(lC) ^ Yi{lD) 



2 



L 



L. 



(3.37) 



11 





Vi{k) Vi{k) 
(B) 



Figure 3: One-loop diagrams for the self-energy of neutrinos in a medium. Diagram (A) contributes only to Vf, 
while diagram (B) contributes equally to Ve, v^i and Vr- 



3.3 Complete vertex 



We have completed the calculation of the diagrams that give rise to the background- ind uced part of the 
proper vertex of neutrinos with the graviton. The complete vertex function defined in Eq. (3.1) i s obta ined 
by ad ding these individual contributions. For the diagrams of Fig. |l| the results are given in Eqs. ( 3.34 ) and 
(3.37), and adding them we obtain 



p/(lA) , w(li3) , p/(lC) , p/(lD) ^ A (W) 



(3.38) 



For the diagrams of Fig. ||, the sum of the individual results given in Eqs. ( 3.13| ), ( |3.27 ) and ( 3.29 ) yields 

(3.39) 



p/(2A) , p/(2B) , p/(2C) , p/(2D) ^ A (Z) 



As we have already mentioned, A^'jT'' contributes to the vertex only when the neutrino is fe, while A'lJl 
contributes equally for all the weak-doublet neutrinos, including the Ve- In addition, the vertex function has 
the tree- level term which, for chiral, left-handed neutrinos is given in Eq. ( 2.18 ). 

In this way we find that the complete effective vertex function for the various neutrino flavors is given by 



iw) 



A 



where 



A„ 



for Ve , 

Notice that the vertex function is symmetric in the indices /x, i^, as it should be. 



(3.40) 



(3.41) 



3.4 Transversality of the vertex 

Before proceeding to the explicit calculations of the neutrino vertex using the one-loop formula given above, 
it is useful to check that the complete effective vertex satisfies the transversality condition 



Here, 



q^'uLik') r(rj(fc, k') UL{k) = UL{k') T^^Jik, k') UL{k) = . 



(3.42) 



(3.43) 
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where lok is the correct dispersion relation for a neutrino mode propagating with momentum K in the 
medium, and ml is the corresponding spinor. As shown several years ago |Q, ^ ^, the spinor ul(^) and the 
corresponding dispersion relation are found by solving the effective Dirac equation 

(|f-I]eff)ML(fc) =0, (3.44) 

where Eeff is the neutrino self-energy in the medium. The chirality of the neutrino interactions dictate that 
Seff has the form 

Eeff = [a]/: + bij)L , (3.45) 

where, in general, a and b are functions of lj}^ and K . In this form, the dispersion relations implied by Eq. 
( 3.44 ) are given by 

UJk = 



1 — a{LUK, K) 



- b{-ujK,K) 

UJK = K - — — , 3.46 

with UJk corresponding to the neutrino and ujk to the antineutrino, which in the general case must be 
considered as implicit equations that must b e solv ed for ujk and ujk as functions of K . In the context of our 
perturbative approach, the solutions to Eq. ( 3.46 ) arc given approximately by 

ujii = K +[l + a{K,K)]h{K,K) 

ujk = K -[l + a{-K,K)]b{-K,K). (3.47) 

At the one-loop level, the neutrino self-energy in the presence of matter is determined by calculating the 
diagrams shown in Fig. which has been carried out in detail in the references cited. The result is that, to 
order the parameter a vanishes while 

be + for Ve 



Vat = \ I ., (3.48) 



where 6e and b'^'^^ are given in Eqs. (3.15) and (3.14), respectively. In this case ujk = K + bmat for the 
neutrinos, and it follows that the spinors satisfy the equation 

(|^-6mat?*)iML(fc) = (3.49) 

and the relation 

fc2 = 2bm^tk -v-b^ (3.50) 

holds. These, together with the analogous relations for u{k') and fc'^, imply the useful formulas 

UL{k') (jiUL{k) = 

e-k'^ = 2b^,tq-v. (3.51) 

We are now in the position to consider the transversality property of the effective vertex. From Eq. ( 3.22| ) 
it follows that 

Q^^lSa = (2p • q - q^)['ip^Pa - "^q^Pa - P^qa + P ■ qVi^a] 

q'^N^^l = {2p-q-q^)ie,o.Xpq^p''. (3.52) 



Then from Eq. (3.35) 



.2 r ^4 



dp 

2^^^^^ = 2M2~''^"'^'' j (2^'^''^^ ~ ml)r]e{p)[2quPa + Puqc - P ■ qVi^a] 

]^be[2qvij + v^i- V ■ q'yu\L , (3.53) 
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and similarly 



'Z'^A^.f ^ = + v4-v q^,]L . (3.54) 



On the other hand, from Eq. (2.19), we obtain 



q^V^^J = + k% + [e - k'')j,]L - -q,[jt + jt']L. (3.55) 



From these, and using t he rel ations given in Eqs. ( 3.49| ) and ( 3.51 ), it is easy to verify that the transversality 



condition stated in Eq. (|3.42| ) is indeed satisfied 



4 Neutrino index of refraction in a gravitational field 

4.1 The self-energy of neutrinos in a static gravitational field 

Our aim in this section is to determine the correction to the neutrino index of refraction in the presence of 
a static gravitational field. To this end, let us consider the scattering of a neutrino by a static gravitational 
potential (/)™*(a;) which is produced by a static mass density p'^^^{x). Defining the Fourier transform of 0°'^* 

by 

with a similar definition for p{q'), the corresponding metric is such that, in momentum space, 

/iA^-(q') = i0(q')(2«V_^M-) ^ (4.2) 

where we have used — 2q^(/) — n^p. The formula in Eq. ( [1.2| ) is the solution to the linearized field equation 
for the metric with the static energy momentum tensor T^"^ ~ v'^v'^p'^^^. Under the influence of such an 
external potential, the off-shell v-v transition amplitude is then 

S,, = -iK{2TT)S{k'> - fc'0)(V;(^)(fc,fc') + A^,(0, Qjjh'^-'ik' - k) . (4.3) 



In Eq. ( [4.3[) we have indicated explicitly the fact that the background part K^^ip., Q) of the vertex function 



does not depend on k and k' separately, but only on the variables Q. and Q that are defined by writing 

= (fc - k'T = (f!, Q) (4.4) 

in the rest frame of the medium. Moreover, we have set = as implied by the delta function in the 
right-hand side. In addition, we now restrict ourselves to a uniform field in space, which means that can 
write 

0(fc' - fc) = (27r)3<5(3) ^j^i _ ;^)^cxt (4 5) 

The justification for this is the usual one. Namely, we assume that we are working in a region of space that 
is microscopically large but macroscopically small, and therefore the external field is approximately constant 
over it. The macroscopic dependence of t he e xternal field on x can then be restored at the end. 
Thus, substituting Eq. (|4.5| ) into Eq. (^), we obtain 

S,, = -i{2nf5^^\k - k')iV^^,\k, k) + A^,(0, Q -> 0)) {2v^^v'' - i^^^) . (4.6) 

Identifying the gravitational contribution to the self-energy by writing 

S,, = -i{2nf5'^^\k - k')llG{k) , (4.7) 

we then find 

Sg = + Sg , (4.8) 
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where 



(4.9) 



As indicated in Eq. (|4.9|) , the gravitational contribution to the self-energy has two parts. We consider first 
Tig{k), which is the contribution from the gravitational field when there is no background medium present. 
Using the expression for V^,y given in Eq. ( |2.19 ), we obtain 



which in terms of the parametrization of the self-energy given in Eq. (3.45) amounts to 



an = 0" 



bg = 0'='''2fc • V . 



(4.10) 



(4.11) 



4.2 The matter-induced gravitational contribution to the self energy 

Since the expression for A^^ contains the factor ^^L^ it is useful to define the vector t\ by writing 

{2v^v^ ~ r^^'^) A^,(r!, Q) = {i^L)t^{n, Q) , 



(4.12) 



in terms of which 

E^(fc) = </)°^VtA(0,Q^O). (4.13) 

In an isotropic background, which we have assumed by writing the distribution functions as in Eq. ( |3.17 ), 
ta can be expressed in the form 



where we have defined 



V = V — 



q- 



(4.14) 
(4.15) 



Since q ■ v — 0, the scalar functions can be calculated by using 

1 



= -^q-t- 

Bg = ^v-t. 

V 

If Ag{0, Q) is not singular in the limit Q ^ 0, then Eqs. ( 1.13| ) and (4.14) imply that 

I]^ = 0-tBG(O,Q^O)^L, 
so that only Eq needs to be evaluated. To prove that this is the case, notice from Eq. ( 3.22 ) that 

{2v^'v'' - r]"") nII\ = q^ {2v''v'' - t]"") N^^J^ = . 
Therefore, Aq and Bq are given by 

Ae+J2fXfAf fori/e, 



(4.16) 



(4.17) 



(4.18) 



Ag = V2Gf 
Bg = V2Gf 



Be+EfXfBf forz^e, 



EfXfBf 



for I^u, I^r 



(4.19) 
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with 



Bf 



(PP 



-2pf -q 
g2 -2pfq 



-(9 



(4.20) 



where we have defined 



h = i'(2v'^v-'~vnNl^Jx- 



From the expression given in Eq. (3.22) for N^^J-^ it is straightforward to see that 

la{0, Q) = (Q' -2P-Q) {8Ej ~ Am) - 2P ■ Q) . 



(4.21) 



(4.22) 



Using this in Eq. ( 4.20 ), it follows that Af{Q, Q) is proportional to the integral of P • Q, which is zero for 
an isotropic distribution. Therefore 

A/(0,Q)=0, (4.23) 

which proves Eq. ( [4.17| ). 

To evaluate i?/(0, Q — > 0), we use Eq. (|^) and the definition of iV^iL given m Eq. (|^) to obtain, 
after rearranging some terms, 



4(0, Q) = 2Ej{8E] - Am} ~ 2Q^) ~ 6Ef{2P -Q-Q^ 



Using this in Eq. ( [4.2QD , 

BfiO, Q) = -i{nj-nj) 

with 



<PP 
(2^ 



F - 2Q\ff - fj) 
-Q2 + 2P- Q 



(Q--Q) 



4{2Ej~m}){ff~ff) 



(4.24) 



(4.25) 



(4.26) 



To evaluate this integral, we make the change of variables P — > P + ^ Q in the first term, and P ^ P — ^ Q 
in the second one, remembering that the distribution functions and F are functions of P. This procedure 
yields 

B/(0,Q^0) = J/ + 0(Q2) , (4.27) 

where 



Jf 



~'S{nf — nj) 



-3(n/ - nj) 



(fP fQ-VpF 
(27r)3 2P- Q 
d^P dF 



{2^Y2Ef dEf 



(4.28) 



In arriving at Eq. ( 1.2g ) we have used the definition in Eq. ( 3.19 ), and in writing the se cond equali ty we 
have used the fact the the function F depend on P only through Ef. From Eqs. ( [1.19| ), ( 4.17 ) and ( 4.27 ) 
we finally obtain 

S^(fc) = (4.29) 



where 



bo 



'V2G 



F X 



Je+J2f^fJf for z^, 



(4.30) 



for l>f^ , Vr 
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4.3 Dispersion relations 



We are now in the position to determine the dispersion relation for the neutrinos in the medium in the 
presence of a static gravitational field. Referring back to Eq. ( 3.45 ), the results we have obtained can be 
summarized by writing 



a 
b 



6mat + bg 



ba- 



(4.31) 



The terms Ug and bg, which are given in Eq. (4.11), represent the contribution from pure gravity and are 
present even in the case that there is no matter. On the other hand, the term fomati given in Eq. ( [3.48 ), arises 
due to the presence of the background medium, independently of whether a gravitational field is present 
or not. Finally, be, given in Eq. (4.30), is the contribution that arises due to the simultaneous presence of 
matter and the external gravitational field. Thus, substituting Eq. (4.31) in Eq. (3.47), and keeping in mind 
that we are allowed to retain only terms that are linear in the external gravitational potential, we obtain 



ujK ^ K + 2K(lf 



(4.32) 



The corresponding formula for ujk is obtained fro m th is by reversing the sign in front of 6mat and bQ. 
Apart from the first term in the right-hand side Eq. ( 4.32| ), which corresponds just to the vacuum dispersion 
relation, the other terms have the following meaning. 

The second term has a purely gravitational origin. If we take it by itself and neglect the effects of the 
medium, it gives rise to the dispersion relation 



(4.33) 



which is equivalent to say that the neutrino and antineutrino acquire an index of refraction given by 

AAo = 1 - 20°'^* . (4.34) 
This is the same result that is obtained by solving the equation 



0, 



(4.35) 



which is the appropriate equation to solve in order to determine the dispersion relation for the photon or a 
massless scalar particle. That purely gravitational term also gives rise to effects analogous to the gravitational 
red shift and bending of light. For example, consider a neutrino beam propagating along the z direction in 
the vicinity of a massive body of mass M situated at the origin, with an impact parameter b. In the eikonal 
approximation, ujKit) and K{x) satisfy the Hamilton equation 



dK 

'di 



— - = -Vujk , 



(4.36) 



which can be used to obtain the transverse component that the momentum develops as a consequence 
of the influence of the gravitational field. Substituting ujk = (1 + '2(j)'^^^)K in Eq. (4.36), we can write 



K 



where the position vector is r = (6, 0, t) and 
for neutrinos 



-2 



dt 



a(/)°'^*(r) 
db 



(4.37) 



i°'^t(r) = —GMjr. From this we obtain the bending formula 
4GM 



K 



(4.38) 



which is the same result as the corresponding formula for photons. 

In contrast to the purely gravitational terms, the matter-dependent terms contribute with opposite signs 
to the neutrino and antineutrino indices of refraction. The term 6niat on the right side of Eq. (4.32), which 
denotes the matter contribution in the absence of a gravitational field, is the usual Wolfenstein term. It 
is flavor dependent, and it is the origin of the MSW mechanism for neutrino oscillations in matter. The 
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remaining terms are the new contributions we set out to determine. One of them is exactly the Wolfenstein 
term multiplied by the external gravitational potential. Since we are considering weak gravitational fields, 
that term is much smaller than the usual Wolfenstein term and can be neglected. Thus, the dispersion 
relation reduces to 

iOK = K + 20^=^* if + 6,nat + ba . (4.39) 
which can be equivalently stated in terms of the index of refraction 

UK 

= (4.40) 

The term that we have denoted by ho is more interesting since it can have a non-trivial dependence on the 
temperature and density of the background material. It is flavor dependent as well, and it arises because 
of the presence of both matter and gravitational fields. In order to be able to consider its possible physical 
effects, in the ensuing section we will estimate its magnitude for various physical conditions of the background 
medium. 

4.4 Estimates of the matter-gravitational effects 

The question we wish to address here is how large are the matter-gravitational effects, represented by the 



term ha in Eq. ( 4.40 ), relative to the Wolfenstein term 6mat- To answer this question we need an estimate 
of the integral if/, which we have carried out in the Appendix, and also make some assumptions about the 
composition of the background for the possible physical situations of interest. 

As an example, let us consider a medium such as the Sun or a supernova. In these , the condi tion s are 
such that the nucleons are non-relativistic. Therefore, using the results given in Eqs. (|A.4| ) and ( [A.8| ), we 
can write for the nucleons [N = p, n) 

-(im^TiN classical nucleon gas 

Jn^{ 3 (4.41) 
degenerate nucleon gas , 



where vfn stands for the Fermi velocity of the nucleon gas. 



Regarding the electrons, if they are relativistic then from Eq. (A. 11) we have Je ~ — Sng. In this case, the 



electron contribution in Eq. ( 4.30| ) is of the order of the Wolfenstein term multiplied by (j)'^^* , and therefore 



it is unimportant. On the other hand, this is not necessarily the case for a non-relativistic electron gas. In 



analogy with Eq. (4.41) for this case 

{—j3mf,ne classical non-relativistic electron gas , 
2^ degenerate non-relativistic electron gas . 

To assess further the possible importance of the matter-gravitational contributions to the index of refrac- 
tion, we also need some knowledge of the magnitude of the gravitational potential that could be involved. 
If we take as a guiding value the potential at the surface of the Sun, 

<^0 = -2 X 10"^ , (4.43) 

we find that the matter-gravitational contributions could be relevant under the appropriate conditions. This 
is particularly true for the nucleon terms, which have the enhancement factors /Jtoat or Vp"^ in the classical 
and the degenerate cases, respectively. This is not unexpected, since the gravitational potential couples more 
strongly to the more massive particles. 



As we have already mentioned, and as Eq. (4.30) clearly indicates, the nucleon contributions to bo are 
the same for all the neutrino flavors, and therefore they are irrelevant for neutrino oscillations involving only 
weak-doublet neutrinos. On the other hand, they are relevant for oscillation phenomena in which sterile 
neutrinos participate and, as we have shown above, they may be important. Furthermore, bmat and be have 
a different dependence on the neutrino coordinate as it propagates through the medium, a property that 
may have also distinctive implications. 
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5 Conclusions 



In this work we have determined the effects of a static gravitational potential on the neutrino index of 
refraction in matter. This has been done by first carrying out the one-loop calculation of the matter-induced 
gravitational couplings of the neutrinos, and then by determining the corrections to the neutrino self-energy 
that such couplings imply in the presence of the gravitational potential. 

As a consistency check of the one-loop formulas, we showed explicitly that the effective gravitational 
vertex of the neutrino is transverse. This required that the correct dispersion relation and wavefunction 
associated with the external neutrino modes be used, and not their vacuum counterpart. We emphasize 
again that it is crucial that, in the calculation of the one-loop diagrams, we have used the full off-shell 
formula for the tree-level gravitational vertex function of the internal fermions in the loop, and not the 
on-shell limit that is customarily quoted. 

As indicated in Section ^, the matter-gravitational contribution to the neutrino index of refraction could 
be relevant in the context of matter-enhanced neutrino oscillations, and in particular in phenomena involving 
the so-called sterile neutrinos. This may occur not only because their magnitude could be non-negligible, but 
also because they have a different dependence on the coordinate compared with the standard Wolfenstein 
term. Whether or not these gravitational effects can lead to interesting observable consequences in specific 
contexts, such as the supernova or the Solar neutrino problem, is an open question that needs further detailed 
study. 

The results presented here indicate that such studies could be well motivated, and our work sets down 
the arena to carry them out on firm grounds and in a systematic fashion. 

Note added: After this paper was submitted for publication, the work by Piriz, Roy and Wudka was 
brought to our attention, in which the tree-level gravitational coupling of the neutrino is also considered 
without violating the Equivalence Principle. However, these authors assume an intrinsic magnetic moment 
of the neutrino, and they concern themselves with the effect of the (vacuum) gravitational interactions on 
the magnetic spin flip oscillations. They do not consider the effect of matter on the gravitational neutrino 
interactions, which is the focus of the present work. 
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A Evaluation of Jf 

We consider the evaluation of the quantity J/ defined in Eq. ( 4.2§| ) for various conditions of the background 
fermion gas. That expression for Jf can be rewritten by taking the derivative of the function F given in Eq. 
( [4.26D and using Eq. ( |3.19| ). In this way Eq. ( [4.28| ) becomes 

Jf^{nf~nj)+Hf (A.l) 

where 



(2^ 



Ef 



f 

2^ 



dE 



f 



■(//-//)■ 



(A.2) 



In order to carry out this integral, we consider the following illustrative cases. 



A.l Classical non-relativistic gas 

In this case, we can put fj — and use 



dEf 



(A.3) 



Using these and remembering that we can approximate Ef ^ mf in the integrand, we then obtain from Eqs. 
( [At|) and ( |a!2|) 

Jf-Hf~ -PrufUf , (A.4) 
where we have used the fact that (3m f ;:!> 1 for a non-relativistic gas. 
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A. 2 Degenerate gas at zero temperature 

For this case, the distribution function is given by 



(A.5) 



where 8 denotes the step function and Ep is the Fermi energy. Thus df /dE = —5{E — Ep). Putting this 
back into Eq. (A. 2) and using Pp = iir^Uf for the Fermi momentum, we obtain 



which in turn imphes 



'f 



1/3 



(A.6) 



(A.7) 



No assumption has been made here about whether or not the gas is non-relativistic. However, in the 
relativistic or non-relativistic limits, this formula reduces to 



Jf 



-5nf relativistic gas 

- non-relativistic gas 



(A., 



where we have used pp = mfvp for the non-relativistic case. 
A. 3 Ultra-relativistic gas 



In this case we neglect the mass of the background particles and therefore we approximate Eq. (A^) by 



Hf =4: 



d^P d 

-PiBiff-ff)- 



(27r)3 dP' 

By carrying out a partial integration and using Eq. (|3.19|) this is equivalent to 

Hf ^ -6(n/ - nj) , 

and therefore 

Jf = -5("/ - nj) ■ 
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